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INFLUENCE  OF  AXIAI.  ENERUY  SPREAD  ON  THE  NEGATIVE-MASS  INSTABILITY 
IN  A RELATIVISTIC  N0NNEUTRA1.  F.- LAYER 


llwan-sup  Uhm 

Depart  meat  of  Physios  and  Astronomy 
University  of  Maryland,  College  Park,  Maryland  20742 

anil 


Ronald  C.  Davidson* 

Division  of  Magnetic  Fusion  Energy 
Energy  Research  and  Development  Administration 
Washington,  D.  C.  20‘>4‘> 


This  paper  investigates  the  influence  of  an  axial  energy  spread  on 
the  negative-mass  instability  in  a relativistic  nonneutral  E-layer  aligned 
parallel  to  a uniform  axial  magnetic  field  , . The  stability  analysis 
is  carried  out  within  the  framework  of  the  linearized  VlaNOV-Maxwel 1 
equations.  it  is  assumed  that  the  E-layer  is  thin  with  radial  thickness 
(2a)  much  smaller  than  the  mean  radius  (R()) , and  that  v/-y^<<l,  where  v 
is  Budker's  parameter  and  y me"  is  the  mean  electron  energy.  Stability 
properties  are  Investigated  for  the  choice  of  electron  distribution 
'function  in  which  all  electrons  have  the  same  value  of  canonical  angular 
momentum  (l'fl-P(1-const . ) and  a step-function  distribution  in  axial  momen- 
tum n . The  neaat tve-mass  urowth  rate  is  calculated  including  the 

1 Z 

important  stablllzine  influence  of  axial  energy  spread  (AE),  and  it  is 

2 ^ 

shown  that  a modest  energy  spread  (AF./y^mc  a tew  pereent  ) is  sulti- 

clent  to  stabilize  perturbat ions  with  axial  wavenumber  satisfying 

'y  i 


*On  leave  of  absence  from  the  University  ot  Maryland,  College  Park.  Md. 
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I . INTRODUCTION 


In  recent  years,  the  properties  of  relativistic  electron  rings  and 


E-layers,  contained  in  a mirror  or  uniform  magnetic  field,  have  been 


studied  experimentally  in  connection  with  plasma  confinement  schemes 


12  3 A 

such  as  Astron  * , single-stage  magnetic  mirrors'5’  , and  electron  ring 


accelerators^  ^ . One  of  the  most  basic  instabilities  that  character- 


izes relativistic  electron  rings  and  E-layers  is  the  negative-mass 


instability  . The  influence  of  a spread  in  canonical  angular  momen- 


tum on  this  instability  has  been  extensively  investigated  in  the  liter- 

Tkn  nvfleAMf-  nanaA  ~ J 1 JL  J - _ , • 


The  present  paper  examines  the  equilibrium  and  negative-mass 


stability  properties  of  a relativistic  nonneutral  E-layer  within  frame- 


work of  the  linearized  Vlasov-Maxwell  equations  , including  the  impor- 


tant influence  of  an  axial  energy  spread.  Recent  experimental  observa- 


tions indicate  that  perturbations  with  finite  axial  wavenumber  are 


significantly  stabilized  by  a small  axial  energy  spread. 


The  present  analysis  is  carried  out  for  an  infinitely  long  E-layer 


aligned  parallel  to  a uniform  magnetic  field  Bn  e [Fig.  1].  Equili- 

U ' vz 


brium  and  stability  properties  are  calculated  for  the  specific  choice 


of  electron  distribution  function  [Eq.  (8)]. 


fe(H’Pe’Pz)  “ 9— 172  6(P0-Po)*(a2-Pz)  *iXmC2r%2 

6 8 2 2irA(2Y0m)1/2  9 0 2 (^mc-U)' 


where  H is  the  energy,  P is  the  canonical  angular  momentum,  p is  the 

0 Z 


— 2 

axial  momentum,  $(x)  is  the  Heaviside  step  function,  U=R+etp  ^-p  2y  is  an 
effective  energy  variable,  and  nQ,  RQ,  A,  yq»  Y and  are  constants. 


Equilibrium  properties  are  examined  in  Section  II.  One  of  the  important 


features  of  the  equilibrium  analysis  is  that  the  equilibrium  distribution 


ii 


i 


function  in  Eq . (8)  corresponds  to  a sharp-boundary  density  profile 
with  uniform  axial  temperature  (Fig.  2).  Stability  properties  are 
investigated  in  Sections  111  and  IV,  Including  the  important  influence 
of  an  axial  energy  spread.  The  stability  analysis  is  fully  electro- 
magnetic and  assumes  that  the  positive  ions  can  be  treated  as  a fixed 
background  ( m ^ ►«•)  on  the  time  scales  of  interest. 

Introducing  the  geometric  factor  g | Eq . (44)]  which  is  a transcen- 
dental function  ot  the  eiRenfrequency  w and  axial  wavenumber  k,  we  obtain 
the  dispersion  relation  iKq.  (5b)] 


(u)-fw  )' 
c 


gw  _ 


(lif-k  Rfl)  - 2k  Rq  - 


A K 


V" 


Here  v is  Budker's  parameter,  AE-A" /2y  m is  the  axial  energy  spread,  4 

is  the  azimuthal  harmonic  number,  to  =eB()/v^mc  is  the  electron  cyclotron 

2 2 

frequency,  and  - 1/y^  is  defined  in  Eq.  (51).  A detailed  numer- 

ical analysis  of  the  dispersion  relation  [Eq.  (5b)  | is  carried  out  in 
Section  IV.  It  is  found  that  the  axial  energy  spread  (AE)  has  a strong 
stabilizing  influence  on  the  negative-mass  instability,  particularly 
when  the  perturbations  have  sufficiently  large  axial  wavenumber  with 
k2R2£l  (Sec.  IV). 

As  a check  on  the  dispersion  relation  given  in  Eq.  (5h>,  it  is 
instructive  to  consider  the  limit  ot  an  ult rare  1 at t v is t ic  (y^>>1), 
inf initisimal lv  thin  (n/R^K))  E-layer,  with  negligibly  small  equilibrium 
self  fields  (p«l)  and  axial  energy  spread  (AE“0).  I n this  case.  Eq.  (5t>) 
can  be  expressed  as 


U-Ju>cr  - -(c2/R2)f(42-k2R2)(v/Y0)gmo/w] 


U.  K.V'Ul  l.l  RRl  I'M  THKORY 


Ah  illustrated  i v>  Kin.  1.  t he  equilibrium  i'o«l  l^ur.it  Ion  consists 
of  a relativistic  nonneutral  K.-laver  that  is  infinite  in  axial  extent 
and  aligned  parallel  to  a uniform  applied  magnetic  tield  Re  < . 
mean  radius  and  radial  thickness  ot  the  K-laver  are  denoted  hv 
and  ,1a,  respectively.  The  mean  mot  ion  ot  the  K- laver  (s  in  t he  azimuthal 
direction,  and  t ho  applied  magnetic  tield  provides  radial  coni  inement 
of  the  electrons.  As  shown  in  Ktg.  1,  cylindrical  point  coot d (nates 
(r,t>,z)  are  introduced.  The  following  are  the  main  assumptions  ^ 
pertaining  to  the  equilibrium  cont igurat ton: 

(al  Kqul  1 Ihr  ivtm  properties  are  independent  ot  z t ' i z 01  and 
azimuthallv  symmetric  ld/A0"01  about  the  c-axls. 

tbl  The  positive  ions  form  an  immobile  tm^  »•>•>,  partially  neut  r.-il  icing 
background.  It  is  assumed  that  t he  equilibrium  ion  density  n'irl 
can  be  expressed  as 


n'(r)-fn*(rl  . I U 

i e 

0 . 

where  t "const,  is  the  t factional  charge  ueut  ra  l l cat  i on . and  n i.  t > 

. • • e 

is  t !*•»«>  1 ect  ron  debs  tty. 

lc)  The  spread  in  radial  and  axial  momentum  t>t  the  electrons  is 
small  in  comparison  with  the  mean  azimuthal  momentum,  i.e.. 


a 


> 0mc 


I 


where  A is  the  chnract er ist ic  spread  in  axial  momentum,  y me' 
azimuthal  electron  energy  at  r*R^,  and  the  liall  thickness  a is 


t."> 


is  the 


related  to  the  spread  in  radial  momentum  bv  Kq.  ilSl. 


(d)  The  electron  motion  is  generally  relativistic,  and  the  mean 
equilibrium  motion  of  the  E-laver  is  in  the  azimuthal  direction,  i.e.. 


d P V " n/r)V,  (>')e0 


where  ( r ) is  the  mean  azimuthal  velocity  of  an  electron  fluid  element 

and  e.  is  a unit  vector  in  the  O-direction. 
vO 

(e)  It  is  further  assumed  that 


1 . 


1 n 


where  v“N  e‘’/mc“  is  Budker's  parameter, 
e 


dr  rn  fri 
e 


(4) 


K 


is  the  number  ot  electrons  per  unit  axial  length  of  the  E-layer, 
c is  the  speed  of  light  in  vacuo , -e  and  m are  the  charge  and  rest 
mass,  respectively,  of  an  electron.  In  F.q.  (4),  and  R,  denote  the 
inner  and  outer  radii,  respectively,  ot  the  E-laver. 

Central  to  a description  of  steady-state  Vlasov  equilibria  are 
the  single-particle  constants  of  the  motion  in  the  equilibrium  field 
configuration.  For  azimuthal ly  symmetric  equilibria  with  3/3z**0, 
there  are  three  single-particle  constants  of  the  motion.  these  are 
the  total  energy  It, 


2 4 1 ">  1 / 2 

tl  = (m  c +c“p‘ ) -e^  . t r) 


thi>  canonical  angular  momentum  1'  , 


*'0  " r|p|rie/c)A0(r)  | . 


i b> 


and  the  axial  canonical  momentum  E 


7 


P - 

z 


P,  . 


(7) 


where  use  has  been  nude  ot  assumption  fd)  . In  F.qs . (5)-(/),  lower 
case  denotes  mechanie.il  momentum,  .f>^(r)  is  the  elect  rosi  at  ic 
potential  tor  the  equilibrium  radial  self-electric  field  |E^(r)» 
-3,*1.,/  ,r|,  and  A^ir)  is  the  0 -component  of  the  equilibrium  vector 


-1. 


potential  |B  ( r ) • r ( S ' > r ) ( rA^l  | . Any  distribution  function  tliit  is 
a 'unction  only  of  tlu  single-particle  onstants  of  the  motion  satistie 
the  steadv-state  Vlasov  |uation.  For  present  purposes,  w<  consider 
tin  equilibrium  distribution  function 

0 n0K0  2 2 

‘>’VV  " — -172  ^Po  Pn)#(A  -p;> 

2r  \ (2YQm) 


18) 


♦(yrnc-U) 
(>  me  “-111  ' 


where  n^  is  electron  dens l tv  at  r»R^ , A , P()  and  t are  constants,  and 
♦Kx)  is  the  Heaviside  st.p  tunc:  ion  defined  by 


| 0 , jcv  0 , 

♦ is)  ■ \ ll)) 

( 1 . X'  0 . 

In  Kq.  (8),  the  en«  rgv  variable  l'  is  defined  bv 

I'-H+tV  (1-p'.  -'>(1m  , HO) 

where  ? ^ is  the  elect t ostat ic  potential  at  the  mean  radius  K,.  i e. . 
i • Note  that  l'  is  a single-particle  constant  of  the  mot  ion 
since  it  is  constructed  iron  a linear  combination  ot  II  and  p*. 

For  a thin  K- layer  equilibrium  consist. nt  with  Fq.  i 'l , 

the  energy  variable  F defined  in  Fq.  (101  can  be  approximated  bv 

i 2 ' 2 

(wit’.  p“<<Y‘m'  c“) 


8 


U:rpr//2YOm+^O^r^+YOmc2  ’ (ID 

where 

i|'0(r)=  [y(r)-Y0]mc‘'-6it.0(r)  (12) 

is  the  "envelope  function",  and 

1/2 

Y(r)mc2={m2cV[c2P0/r+(e/c)A0(r)  ]2}  , (13) 

is  the  electron  energy  associated  with  the  azimuthal  motion.  Moreover, 

2 2 

6 4>q ( r ) is  defined  by  6<{>Q(r)=<}ip(r)-<()Q.  and  YqIhc  =Y(r=RQ)mc  is  the  azimuthal 
electron  energy  at  r=Rg.  The  envelope  function  ^(r)  can  be  further 
simplified  for  a thin  E-layer  satisfying  Eqs.  (2)  and  (3)  by  Taylor 
expanding  Eq.  (12)  about  the  mean  radius  Rq . Retaining  terms  to 
quadratic  order  in 


P 


=r-R„ 


(14) 


and  making  use  of  equilibrium  radial  force  balance  ’ on  an  electron 
fluid  element  at  r=Rp,  we  find 


Vr)  = 1 Y0m‘,rp2 


(15) 


where 


2 2 2 ,1  2,.  . 2 
Wr=Wc4Wp(1_Y0f)/Y0 


(16) 


2 0 2 

is  the  radial  betatron  frequency,  w =47rn  (Rn)e  /ynm  is  the  plasma 

p e U U 

frequency-squared  at  >"=Rq,  and  ioc=eB^/YQmc  is  the  electron  cyclotron 
frequency.  [See  Refs.  16  and  19  for  a detailed  discussion  of  Eq.  (15)]. 
Making  use  of  Eqs.  (11)  and  (15),  several  interesting  properties 


can  be  deduced  for  the  class  of  thin  E- layer  equilibria  described  by 
Eq.  (8).  For  example,  it  is  straightforward  to  show  that  the  electron 


9 


density  profile  can  be  expressed  as 


0 3 n R0  2 2 

n^(r)=  d pfe(H,P  ,Pz)=n0  — *P  ) . 


where 


,-r  ,2,  2,1/2 

a=  [ 2 ( y — y q)  c /YqCJjJ 


(H) 


(18) 


is  the  half-thickness  of  the  E-layer.  In  obtaining  Eq.  (17),  use 
has  been  made  of  (3/3p  )P  =r  and  the  identity 

0 o 

dx(a2-x2)"1/2=7r'2  . (19) 

. 

0 

Since  the  E-layer  is  assumed  to  be  thin,  we  approximate  Rq/t=1  in 
Eq.  (17),  and  the  electron  density  profile  reduces  to 


0,  ^ 
n (r)  = 
e 


*0  • Ir-R0l<a  ’ 

0 , otherwise  . 

Moreover,  it  is  straightforward  to  show  that  the  axial  electron 
temperature  profile  can  be  expressed  as 


(20) 


T°(r)= 


( ,3  , 2,  ,f0 

d p(Pz/Y0m)f 


A 


,3  ,0  . 2 

d pf  = A / JY„m=con.-it . 
e 0 


(21) 


for  |r-Rg|<a.  The  electron  density  and  axial  temperature  profiles  are 
illustrated  in  Fig.  2.  Evidently,  the  equilibrium-  distribution 
function  in  Eq.  (8)  corresponds  to  a sharp-boundary  density  profile 
[Eq.  (20)]  with  uniform  axial  temperature  [Eq.  (21)]. 

Finallv , we  conclude  this  section  by  noting  that  the  spread  in 
axial  momentum  in  Eq.  (8)  yields  a corresponding  spread  in  total  energy 
H.  From  Eqs.  (10),  (11).  and  (15),  we  express 


Y0mc  - e 


to  + 2^»  + 


2 2 


-x + n ynmi>  p 

2Yr.ni  2 0 r 


(22) 


Within  the  context  of  F.q.  (2).  it  is  valid  to  approximate  Kq . (22)  bv 


H^0mcJ-ei0+P;/2Y0m  . 

It  is  evident  from  Kq.  (8)  that  2 A is  a measure  of  the  total  spread 
In  axial  momentum.  Therefore,  from  Kq.  (23),  the  total  energy  spread 


2 

AE-  A *72vrtm 


(2 


can  be  approximated  by 


11 


where 


^(x)=-e  d p v f (x,J>) 


is  the  perturbed  current  density. 


r° 

= ej 

J — a 


dTexp{-iu)x}  \E(x’)  4-  — 


?'*«<*’>)  a o 


ap'  e 


is  the  perturbed  distribution  function,  and  r-t'-t.  To  make  the 
theoretical  analysis  tractable,  we  Fourier  decompose  the  0-  and  in- 
dependence of  all  perturbed  quantities  according  to 


ij <(x)=I  dki  (k,r)exp{  i ( £0+kz)  } , 

^ oJ  1 


III.  ELECTROMAGNETIC  STABILITY  ANALYSIS 


A.  Linearized  Vlasov-Maxwell  Equations 


In  this  section,  we  examine  the  linearized  Vlasov-Maxwell 
equations  for  perturbations  about  a thin  E-layer  equilibrium  described 
by  Eq.  (8).  To  determine  stability  properties,  we  adopt  a normal- 
mode approach  in  which  all  perturbations  are  assumed  to  vary  with  time 
according  to 


6iK^,t)=^(^)exp{-iut}  , 


with  lmw>0.  Then  the  Maxwell  equations  for  the  perturbed  electric 
and  magnetic  fields  become^ 


yxKx) 

'V  % 'Xj 


VxB(x) 

% , <\j 


i - B (x)  , 

C *1*  *\j 


f - w - 1 7 w ’ 


where  k is  the  axial  wavenumber  and  1 is  an  integer.  Making  use  ol 


Ki|s.  (2r>)  and  (28),  II  Is  st  ra  1 ^ht  forward  t o show  (hat 


h , (k  , r ) - 1 V'  p* 
r f z rut 


*(r)  iVl  Kfz(k,r) 
P 


will*  VC 


**  ) > ') 

*-  *•  / * * *■ 
' -u»  fc  -k 


and  tho  function  4*  (.  r f is  deft  nod  hy 


if  i r)-lrE  (k,r)/ i . 


Equation  (29)  can  bo  expressed  as 


A .(k.r)  - l'  lif(r)  * ~ E tk.r) 

• MV.  Ilil  ( . t . 

P 


Ik 


outside  tho  E-layer,  because  tho  perturbed  azimuthal  current  density 


d (k,r)“0  in  the  vacuum  region. 


The  perturbed  axial  and  azimuthal  electric  Holds  |K,ik,r)  and 


E (k,r)l  are  continuous  across  the  beam  boundaries  ir  Kj  and  t -K,), 


as  Is  the  function  if  ( r ) (Fq.  (U)|.  Integrating  Eq . f.’1')  t row  t -U 


to  r“K  ,+iS  and  taking  the  limit  *0^  , we  ohta 


I u 


B ( k , R.,)  U y k , K , > 

< z 2 < / I 


\ in')  if  ( K ) 

i It 


c v’ 

♦ I' 

10 


ik 


if  ( r ) / r * , F , ik ,r ) 

P*« 


lb 


where  the  et  feet  Ive  susceptibility  is  del  tneil  bv 

. t 


\ in') 'f  I K j) 


dr  d (t  ) , 

t 0 


i 


and  if(R  ) dettotes  lint  ( K *iS).  For  presettt  purposes,  it  is  alsi 


assumed  that 


I m—  C.mt,  | <<o)r  t 


t(ui  /hi  )<--K.,/a  . 
c r 0 


Within  the  context  of  assumption  (d)  and  Kq . ( 33) . it  is  valid  to 

. _ 14,1b 

approximate 


i(r)^|(K0)  . 


Vk-r)'E.,<k-V  • 


",o(k'r’ -"<o<k'V 


For  convenience,  in  the  subsequent  analysis  we  Introduce  the 
normalized  electric  and  magnetic  wave  admittances,'^*^1  d and  b , 

ii 

defined  at  the  inner  and  outer  surfaces  of  the  K- layer  hv 


d+— |r(3/nr>Ksjs(k.r)  |R+^|  (E^lk.K,)  1 


d_*  | r (3  /3  r)  F.(^  (k,  r)  1R_/|  < f.(?(k.K)) 


b+“  - ( R ^ ( k , K ^ ) / l r ( 3 / 3 r ) R ^ _ ( k , r)  1R 


b *(B,  (k,R.)/(r(3/3r)Il  (k,r)|u 

v 7-  1 v ~ K I 

For  a beam  in  vacuum,  the  values  of  b(  and  d depend  on  the  geometric 
configuration  (Appendix  ).  Making  use  oi  Kqs.  (V)-(lb),  it  is 
straightforward  to  show  that 


(b_+h+)  |*  (R0)+ikij?,(k.K0)/p‘  1 - v(,.i)^(K0) 


14 


Fourier-decomposing  the  perturbed  fields  in  Eq.  (25),  we  obtain 


B£Q(k,r)  2 B^Oc.r)  + 2 3r  ^(k.r) 

P r cp 


(40) 


in  the  vacuum  region.  Since  B^(k,r)  is  continuous  across  the  beam 


boundaries  and  approximately  constant  within  the  beam  (Eq.  (36)], 
it  is  straightforwrad  to  show  from  Eq.  (40)  that 


1 kf^  Ur(a/3r)Eez(k,r)]  +-[r(3/3r)E^(k,r)  ] _} 


(41) 


=Bez(k«R2)-Bez(k*Ri)  • 


Substituting  Eqs.  (33)  and  (37)  into  Eq.  (41),  and  making  use  of  Eq. 
( 35) , we  f ind 


4nk 


(d_+d+)Eez(k,R0)  = -i  V1  X <<*>)♦  (ty 


(42) 


After  some  simple  algebraic  manipulation  that  makes  use  of  Eqs.  (39) 


and  (42),  the  function  <£(Rq)  can  be  expressed  as 


2 li 


gx(w)  = 1 , 


(43) 


where  the  geometric  factor  g is  defined  by 

,2 


R 2 


2,  2 

w /c 


b +b+ 


d_+d+ 


(44) 


Evidently,  an  evaluation  of  the  effective  susceptibility  x(m)  *s 
required  for  a detailed  stability  analysis. 


B.  Ef feet ive  Suscept lbl 1 lty 


In  this  section,  we  evaluate  the  perturbed  azimuthal  current 
density,  and  the  effective  susceptibility  x(t»)  defined  in  Eq.  (34), 


I'l 


From  the  Maxwell  equation  ^*£(^)  i (u>/c)  . wo  obtain 

B|JS(k.r)  - J-“  «er(k,r)-ikf(r)  . 

E (k,r)  + r B (k  ,r)  - - 1 Mr)  . 

V r lo  (z  i r y 


t45) 


Since  t ho  elgent roquoiu-v  tor  the  negnt  i ve-mass  instability  l'>  given 

approx  imatelv  hv  i,>  <<o  lor  V*\r)«0,  and  sinoe  the  thickness  of  K 

c z 

layer  Is  much  smaller  than  its  major  radius,  Eq.  (45)  can  be  approximate 
by 


E^(k,r)-(v  /i  lB  ^ (k , r >»- i k$ ( r)  , 


E ^ (k , i >+(v.)  \ i B >(k,r)*-(.>/;»r)^(r) 


( .M 


After  some  straightforward  algebra  that  makes  usi  of  Kqs.  (2/)  and 
(46),  we  obtain 


0 


fci(x.^)— e dT  exp{-i(o i } V ' ^ (x ' ) • f^H.P  ,PZ1  , (47) 

—OP  ' f\, 


where  $(x)“$  (r)  exp{  1(  f.0+kz)  } , and  leading-order  terms  have  been 
retained  [see  Fq . (15)].  Fourier  decomposing  Eq.  (47)  according  to 
Eq.  (281  and  making  use  of  the  appioxtmatlon  j ( r)  2 (K())  [Eq.  v '<>)  ) < 
the  perturbed  distribution  function  can  be  expressed  as 


fp,(r.^)  — 1e4,(R0)I|k(.l/liP  UcO/aP  ,!  ir(H,P  P?) 


(4:.) 


where  the  orbit  integral  I i?  defined  bv 


,-n 


1“  di  expl i [ t (o ’-0) I k(z ' -z) -a 


(441 


As  indicated  in  Eq.  (49),  the  particle  trajectories,  0 ' ( i ) and 


lb 


z'(t),  In  the  equilibrium  fields  are  required  In  order  to  evaluate 
the  orbit  Integral  I.  Making  use  of  the  Hamiltonian  in  Kq.  (5), 

It  Is  straightforward  to  obtain  the  equations  of  motion.  Assuming 
that  the  electron  orbit  passes  through  the  phase-space  point  (z,p  ) 
and  0'»b  at  time  t'"t,  we  find 


z'- z+p2r/y0m  , 


O'-O+Cco^-pdiP^/Y^inR^)  r . 


(SO) 


where 


/ui  - 1 / > 


0 * 


(51) 


and  M’ (»1’  -1'^.  [Set*  Ket  . 1 (>  tor  a detailed  derivation  of  Kq.  (50)1. 
Substituting  Kq . 1 50)  into  Kq.  (4l>)  , we  obtain 


1 - i ( - -kp^/ ( / -y^mR^) 


2,-1 


(52) 


From  Kq.  (2b>,  the  perturbed  azimuthal  current  density  can  be 
expressed  as 


f 1 ' 

.1  (k.r)--ew  Rn  d p t ,(r,n) 

i 0 c 0 e { 5 


(53) 


where  we  have  uppt oxlmated  v R^.c  for  a thin,  relativistic  K-laver. 
Substituting  Kqs.  (48)  and  (42)  into  Kq.  (53),  and  carrying  out  the 
momentum  integration,  the  perturbed  azimuthal  current  density  can  be 

expressed  as 


> 3 > > -3 


; . s x "Vc  ‘,r-k‘Ro 

■'  tk  .t  > - -if  (K  .)  , ,, 

'o  0 (ui -tw  ) -k  A />,,m 
c U 

where  use  has  been  made  of  a K^.  The  number  of  electrons  per  unit 
>x'  '1  1 earth  ot  the  1 -lavet  (N  f can  be  determined  bv  substituting. 


t 4a  f 
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Eq.  (17)  into  Eq.  (4).  Within  the  context  of  Eq.  (2),  we  obtain  the 
approximate  result  N^^itan^RQ.  Eliminating  nQ  in  favor  of  N in 

Eq . (54)  and  substituting  Eq . (54)  into  Eq.  (34)  gives 

2 9 2 2 

Nee  wc  -k“Rn 

X(w)  j 2 — 2 — 2 — TT  • 

2iry0mR0  (a >-Iuq)  -lc  A /yV 

Equation  (55)  is  used  in  Sec. IV  to  complete  the  stability  analysis. 


(55) 


IV.  NEGATIVE-MASS  STABILITY  PROPERTIES 


Substituting  Eq.  (55)  into  Eq.  (43)  and  making  use  of  Eq . (24) 

2 2 

and  the  definition  of  Budker's  parameter,  v-N^e  /me  , the  dispersion 
relation  for  the  negative-mass  instability  can  be  expressed  as 


, , v 2 C2  V gwc  , .2  ,2  2.  2„2  AE 

(“"*u,c)  = - ^ y:  — (lie  “k  V ’ 2k  Ro  — 2 
R0  K 0 Yomc 


for  l<fc<Rg/a.  In  Eq.  (56),  g is  the  geometric  factor  defined  in 

Eq.  (44),  AE“  A /ly^m  is  the  energy  spread  [Eq.  (21)],  and  -l/y^ 

[Eq.  (51)].  It  should  be  noted  that  Eq.  (56)  is  valid  only  when 

|u>-lu>c|<<  ii>r  [Eq.  (35)].  The  growth  rate  (D^-Imu  and  real  oscillation 

frequency  $1  *Rew  can  be  determined  from  Eq . (56)  for  a broad  range  of 

2 

I»rameters  k,  v/y^  and  a-AE/y^mc  , by  solving  numerically  the  full  trans- 
cendental dispersion  relation.  [Note  that  the  geometric  factor  g in 
Eq.  (56)  is  generally  a complicated  function  of  the  complex  eigen- 
f requency  w. ] 

For  present  purposes,  to  illustrate  the  influence  of  the  axial 
energy  spread  on  the  negative-mass  instability,  we  calculate  the  normal- 
ized electric  and  magnetic  wave  admittances  defined  in  Eqs.  (A. 3)  and 
(A. 4)  by  approximating  in  the  expressions  for  d(  and  bt.  Figure 

3 is  a plot  of  sum  of  the  wave  admittances  (b_+b+)  and  (d_+d+)  versus 
normalized  axial  wavenumber  kR  , for  f=6,  yr,=5,  a/R  -0.02,  T_/R  *1.5, 


and  T^=0  (no  inner  conductor) . In  Fig.  3,  the  solid  and  broken  curves 
correspond  to  (b  +b+)  and  (d  +d  ) , respectively.  Evidently  (b_+b+)  and 
(d  +d+)  are  slowly  varying  functions  of  kR^  for  the  range  OikR^^.2,  where 
the  approximation  is  certainly  valid.  Negative-mass  stability 

properties  are  illustrated  in  Figs.  4 and  5 for  the  range  of  axial  wave- 
numbers  0£kRQ^,2.  Generally  speaking,  the  quantities  (b_+b+)  and  (d_+d+) 
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either  vanish  or  become  Infinite  at  certain  values  of  k.  [For  example, 

from  Fig.  3,  (b_+b+)-0  at  kR0=3.1,  and  (b_+b+)-*®  at  kRQ=3.5.] 

Normal  inodes  with  axial  wavenumber  k for  which  (b_+b+)-*4)  or  (d  +d+)-»-0, 

are  purely  transverse  electric  (TE)  or  transverse  magnetic  (TM) 

waveguide  modes,  respectively.  The  dispersion  relation  [Eq.  (56)] 

18 

for  the  purelv  TE  and  TM  modes  will  be  investigated  elsewhere. 

Figure  4 shows  a plot  of  normalized  growth  rate  ui./ui  versus 

i c 

kR0  obtained  from  Eq.  (56)  for  f-0.1,  v/yq-0.02,  several  values  of 
a-AE/-*0mc',  and  parameters  otherwise  similar  to  Fig.  j.  Two  important 
features  are  evident  from  Fig.  4.  First,  the  axial  energy  spread 
stabilizes  perturbations  with  axial  wavenumber  larger  than  a certain 
critical  value.  For  example,  perturbations  with  kRg^l.03  are  completely 
stabilized  by  a 4t  energy  spread  (a”0.04  in  Fig.  4).  Second,  for  given 
kR^»*0,  the  growth  rate  can  be  substantially  reduced  hv  increasing 
the  energy  spread.  Those  results  are  consistent  with  recent  experimental 
observations  bv  Destler  et  al . , who  find  that  a finite-length  E- 
laver  (in  which  kK(1  is  necessarily  nonzero)  is  stable  whenever  the 
energy  spread  is  increased  by  a sufficiently  large  amount. 

Stability  boundaries  in  the  parameter  space  (kR^,  AE/>^mc*’) 
are  illustrated  in  Fig.  5.  The  solid  curves  correspond  to  the  stability 
boundaries  (<.>^”0)  obtained  from  F.q.  (5b)  for  several  values  of  v/> ^ 
and  parameters  otherwise  similar  to  Figs.  3 and  4.  For  a given  value 
of  v/> (1 , the  region  of  lkR(1,  AF./>  ^me")  parameter  space  above  the  curve 
corresponds  to  stability  (u>j«0),  whereas  the  region  of  parameter 
space  below  the  iirve  corresponds  to  instability  (u'j-O).  For  sufficiently 
low  beam  density  (u/>^«0.01,  sav) , it  is  evident  from  Fig.  5 that 
the  system  is  stable  provided  kR^'- 1 and  a^O.02.  However,  even  for  low 
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beam  densities.  It  is  evident  from  Fig.  5 that  long  wavelength 
perturbations  with  kRQ<0.5  are  not  stabilized  by  an  axial  energy 
spread,  at  least  for  the  choice  of  distribution  function  in  Eq.  (8). 

In  this  regard,  we  hasten  to  point  out  that  high-harmonic  pertur- 
bations with  small  axial  wavenumber  can  be  stabilized  by  a spread  in 
canonical  angular  momentum.^ 

We  conclude  this  section  by  noting  that  stability  properties  also 
exhibit  a sensitive  dependence  on  the  location  of  the  conducting  walls. 
Moreover,  the  equilibrium  and  stability  analysis  presented  here  can 
be  extended  in  a straightforward  manner  to  hollow  electron  beams 
that  propagate  along  the  z-axis  with  a nonzero  axial  velocity  . 


18 


- 1 

f* 
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V . S UMMARY  ANT)  CON  Cl . L'SJ  ON  S 

In  this  paper,  we  have  examined  the  inlluence  ol  axial  energy 

spread  on  the  negative-mass  install!  lity  In  a relativist ie  nonneutral 

E-layer.  The  analysis  was  earried  out  within  the  I ramework  of  a 

linearized  VI asov-Maxwel 1 equation.  The  equilibrium  (See.  ID 

and  negative-mass  stability  (Sees,  til  and  IV)  properties  were 

investigated  in  detail  tor  the  choice  of  distribution  function 

in  which  all  electrons  have  the  same  value  of  canonical  angular  momentum 

and  a step-function  distribution  in  axial  momentum  p [Eq.  (8)]. 

A detailed  numerical  analysis  of  the  dispersion  relation  [Kq.  (5b)  | 

was  presented  in  Sec.  )Y.  One  of  the  most  important  conclusions 

ot  this  study  is  that  an  axial  energy  spread  can  have  a large  influence 

on  stability  behavior.  in  part icular,  perturbations  with  suf f icient ly 

large  axial  wavenumber  (kR^l)  can  be  completely  stabilized  by  a small 

axial  energy  spread  AbtfO.  Moreover,  in  the  special  limiting  ease  when 

A IX)  and  icl.  the  stability  properties  are  consistent  with  the  results 

1 A 

previously  obtained  bv  Briggs  and  Neil. 

Finally,  we  emphasize  that  the  numerical  analysis  presented  in 
Sec.  IV  is  based  on  the  assumption  that  the  geometric  factor  g is  a 
slow  varying  function  ot  k.  Although  this  is  a reasonable  approximat ion 
for  ()■  k!Cj,<2  (Fig.  '].  we  expect  significant  mod  if  leaf  ions  to  the  stability 
behavior  whenever  (h  4b+)  v0  or  (d  +d()v0.  The  excitation  of  transverse 
electric  lib  +b('M)]  and  transverse  magnetic  [id  +d+)‘0]  waveguide 
modes  bv  the  negat t ve- mass  instability  is  currently  under  invest igat ion . * ' 
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APPENl'1  X 
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WAVE  ADMITTANCE  AT  THE  BOl’NDARlES  OF  AN  K- LAYER  IN  A 
CYLINDRICAL  CAVITY 


In  this  section,  we  obtain  expressions  for  the  wave  admittances 
at  the  boundaries  ot  an  E-layer  in  a cylindrical  cavity.  Since  the 
perturbed  current  density  vanishes  in  the  vacuum  regions  outside  the 
E-laver,  the  Maxwell  equations  in  this  region  can  be  expressed  as 


/la  a t2  2 \ 

l r dr  r u - 2 ^ P > 


E,  (k.r) 


S,  (k.rl 
(z 


(A.n 


In  obtaining  Eq.  (A.  1) , use  lias  been  made  of  Eqs.  (25),  (.28)  and  i H>)  . 
The  solution  tv'  Eq.  (A.l)  is  given  by 


1 E (k  , r) 

1 1 

[ h i k , r ) 
(z 

rl 

Ad  (pr)+HN  tprl  , 

< V 

Cd  ,(pr)+DN  fprl  , 
i C 


R rvT 


0 ’ 


Vr.Rt 


( A . 2 ) 


where  .1  ,(pr)  and  N(.  iprl  are  Bessel  functions  of  the  first  and  second 

kind,  respectively,  and  T ^ and  are  the  radii  v’t  the  inner  and 

outer  conductors , respectively  (see  Eig.  1).  In  Eq.  (A. 2),  A,  B, 

0 anil  D are  constants  which  are  determined  from  Eqs.  (.17)  and  t,  18 1 

and  the  boundary  conditions  for  the  perturbed  fields. 

Making  use  ot  the  boundary  conditions  E,  tk,T  l-E,  (k,l‘  )=0, 

vz  t vz  0 

t tie  electric  wave  admittances  at  the  inner  and  outer  boundaries  ot 
E-laver  [Eq.  can  bo  expressed  as 

K,P  d’.lpR^N^pT^-d  ,(pT0)N'c(pR:) 

d+  < •*  e(pR  ,)N  ((pT0)-d  t(pT0^N  e(pR2) 

VA.  .O 

RjP  d’ipR^N  .tpT^-d  ^pT^N’ipK,! 
d-  ” i d APK^N  Ap'Y-d  ApTjVn  c\pK,  l * 
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FIGURE  CAPTIONS 


Fig.  1 Equilibrium  configuration  and  coordinate  system. 

FIr.  2 Electron  density  |F.q.  (20)]  and  axial  temperature  profiles 
[Eq.  (21)]. 

FIr.  3 Plot  of  sum  of  wave  admittances  (b  +b+)  and  (d  +d+) 

(Appendix)  versus  normalized  axial  wavenumber  kR^ , for 
<-6,  Yq=5,  a/R^*=0.02,  T(1/R(1=1.5  and  Tj»0. 

FIr.  4 Plot  iif  normalized  growth  rate  versus  kK^  [Eq.  (5b)], 

for  f-0.1,  v/y  =0.02,  and  several  values  of  a*AE/y_mc‘~ . 
Parameters  are  otherwise  similar  to  Fig.  3. 

Fig.  5 Stability  boundaries  [Eq.  (5b) [ in  the  parameter  space 

(kR(1>  AF./y^mc* ) for  several  values  of  v/y^  and  parameters 
otherwise  similar  to  Figs.  \ and  4. 


